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The energy and momentum dependence of the gluon self-energy is investigated in a color 
superconductor with two flavors of massless quarks. The presence of a color-superconducting 
quark-quark condensate modifies the gluon self-energy for energies which are of the order of 
the gap parameter. For gluon energies much larger than the gap, the self-energy assumes 
the form given by the standard hard-dense loop approximation. It is shown that this mo- 
dification of the gluon self-energy does not affect the magnitude of the gap to leading and 
subleading order in the weak-coupling limit. 

I. INTRODUCTION 

Single-gluon exchange between two quarks is attractive in the color-antitriplet channel. Therefore, suf- 
ficiently cold and dense quark matter is a color superconductor 0|. When the quark-chemical potential 
H ^ Aqcd, asymptotic freedom |^ implies that the strong coupling constant g at the scale /i is much 
smaller than unity, g{p) «C 1. This allows a controlled calculation of the color-superconducting gap param- 
eter (j) in the weak-coupling limit. 

In order to compute the gap parameter, one has to solve a gap equation. In general, the gap equation 
determines the gap parameter (j){K) as a function of 4- momentum if^ = (fco,k) of the fcrmionic quasipar- 
ticles. However, of physical interest is only the value of the gap function on the quasiparticle mass shell 
[H, i.e., for fco = Cfc, where ek = y/ (ijJ^ — ^i)^ + 0^ is the fermionic quasiparticle excitation energy. Here, uj^ 
is the kinetic energy of non-interacting particles, uj'^ — k = |k| in the ultrarelativistic case. Then, at zero 
temperature, T = 0, the gap equation is typically of the form 

(/)(efc,k) JyJ dcos6'A(efe-eg,k-q),^(e,,q) . (1) 

Here, p{^) is the density of states at the Fermi surface, ^ /i^ in the ultrarelativistic case. Furthermore, 
it is assumed that the attractive interaction between fermions is mediated by boson exchange, with g being 
the fermion-boson coupling constant and A(P) the boson propagator; 6 = arccos(k, q) is the angle between 
the 3-monienta of in- and out-going fermions in a scattering process. 

For an ordinary superconductor, standard BCS theory generally assumes that the attractive phonon- 
exchange interaction between electrons is local, i.e., a point-like four-fermion interaction, and isotropic, 
A(P) = 1/A^, where A is a constant with the dimension of energy. As a consequence, the integral over cos6' 
in Eq. (|^) is trivial. In a superconductor, the maximum phonon momentum is po, the Debye momentum. 
In weak coupling, pu gfi <^ fi. Consequently, the momentum exchange between electrons in a phonon- 
mediated scattering process is limited to the range |k — q| < po- Let us consider k = kp, where kp is, the 
Fermi momentum, and denote the value of the gap function at the Fermi surface by (p. The restriction on 
possible phonon momenta then translates to limiting the g-integration to the (in weak coupling, narrow) 
region kp ~ Pd i£ Q i£ kp + po around the Fermi surface. Assuming that the momentum dependence of 
0(eg,q) near the Fermi surface, q ~ kp, is weak, (/)(eg,q) (j) — const., the gap equation simplifies to 

(f, = In ( + + '^^ 

CBcs y 
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where 1/cbcs ~ PilA/^^- Here, lud = vfPd is the Debye frequency and vp the Fermi velocity, ■yp = 1 in 
the ultrarelativistic hmit. For any, even arbitrarily weak, attractive interaction there is always a nontrivial 
solution ^ ^ to this equation. Consequently, the constant factor common to both sides can be divided out. 
In weak coupling, 1, and in order to solve Eq. (H) the logarithm has to be ~ l/g^ ^ 1, to compensate 

the small prefactor <C 1. This requires (j) <C wd- Then one may approximate ln[(ci;£) + y'uj^y + cj)^) / (j)] ~ 
ln[2uj£i/(j)), and the solution is 

= &BCS M exp ^-^^1^^ , (3) 

where fescs = Swd//^ is a dimensionless constant. This discussion shows that the constant cbcs in the 
exponent is determined by the constant in front of the "BCS logarithm" J dq/eq ln(6BCS l^/4')j while the 
constant in the prefactor, 6bcSi is determined by the constant "under" the BCS logarithm, in this case the 
size of the g-integration region. 

The discussion can be easily generalized to non-local interactions of finite range, for instance massive 
scalar boson exchange Q|. In this case, A(P) = {Mg — P^)~^ is the boson propagator, with Mb the boson 
mass, = Pq — p^. For the following, let us assume that fermions are massless, and that the boson mass is 
generated by in- medium effects, Mb — i.e., in weak coupling the different energy scales in the problem 
are ordered, 4> ^ The angular integral is no longer trivial; in the ultrarelativistic case and near 

the Fermi surface, k q kp = 

^(^)^(^)^ 2(l-cos.)+Ml/M- ' 

where the energy dependence of the boson propagator has been neglected, since po = Ck — ^ (j) Mb- 
Obviously, the term enhances the contribution from small-angle fermion-fermion scattering, cos 9 ~ 1, to 
the gap equation (|l]). Performing the angular integral introduces an additional logarithm, in the following 
called "collinear logarithm" , in comparison to Eq. , 

where I exploited the hierarchy of scales in weak coupling. The solution is rather similar to the one in BCS 
theory, Eq. (|), 

Quantitatively, the difference to Eq. (^) is that, in weak coupling, the gap becomes larger, since the coupling 
constant is replaced by g'^\n{2 h/Mb) ^ g'^\n{l/g) » g"^ . 

In QCD, gluon exchange is also a non-local interaction. The difference to the previous case is that gluons 
are massless and thus the interaction has infinite range. In the vacuum, this holds for electric as well as 
magnetic gluons. In a dense medium, however, electric and non-static magnetic gluon exchange is screened. 
Essentially, these gluons behave like the scalar bosons of the previous example with Mb = mg, where 



Nf 

is the gluon mass at T = 0, with Nf being the number of massless quark flavors participating in screening 
electric and non-static magnetic gluons. Screening reduces the range of the interaction to distances < iti~^. 
On the other hand, the exchange of almost static magnetic gluons is still unscreened and thus of infinite 
range. 

In the vacuum, the gluon propagator is A(P) ^ — 1/P^, and the previously discussed collinear enhance- 
ment in the gap equation is much stronger: instead of being cut off by a mass term Mb, it is cut off by 
the gluon energy, po = Cfc — Eg ~ Mb — gfJ-- (In the previous case, the boson energy was negligible 
compared to its mass.) The collinear logarithm in the gap equation is now ~ ln[2 /i/|efc — e^j] ~ ln(2 fJ-/(t>). It 
is thus of the same order as the BCS logarithm, and consequently much larger than the collinear logarithm 
\n{2 (x/Mb) ~ ln(l/f;) in the previous case. In a dense medium, the gluon propagator is more complicated. 
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but since the exchange of almost static magnetic ghions is unscreened, it still gives rise to the same collinear 
logarithm ln(2 n/ 4>) as in the vacuum 

In QCD, there is no restriction on the magnitude of the gluon momentum in the gap equation. Nevertheless, 
in weak coupling the gap function (/)(eg, q) is strongly peaked at the Fermi surface q = kp Effectively, 
this again restricts the q-integration to a narrow range kp — S < q < kp + S around the Fermi surface, 
(f) ^ S ^ fi. For an order-of-magnitude estimate (but not for a numerically correct solution), in this range 
one may take (/>(efc, k) ~ 4>{£q, q) ~ = const., such that the gap equation assumes the form 



9' 



(8) 



The first term contains two powers of the logarithm \n{2 6/(j)). The first one is the BCS logarithm. The 
second one is the collinear logarithm arising from the exchange of massless gluons. (Any constant under 
this logarithm different from 2S has been absorbed in b'.) The value of the constant is determined by 
how many gluon degrees of freedom are unscreened and thus give rise to a collinear logarithm in addition 
to the BCS logarithm. In vacuum, these are all gluons, while in medium, only almost static magnetic 
gluons contribute to c^. The second term contains only one power of the logarithm ln(2 5/0): the BCS 
logarithm. In a strongly interacting, dense medium, this logarithm arises from exchange of screened electric 
and non-static magnetic gluons The value of b' is determined by the number of degrees of freedom of 
such gluons. As seen above, screened boson exchange also gives rise to a collinear logarithm, but this time 
of the form ln(2 h/Mb), Mb = mg for gluons. This logarithm has been absorbed in the constant b' . Finally, 
all contributions without any logarithm have been summed up into the constant d. 

In weak coupling, <C (5 and hi{2S/(j)) 3> 1. Thus, all unscreened gluons contribute to leading order to 
the gap equation, i.e., to the first term in Eq. (^, while all screened gluons contribute to subleading order, 
i.e., to the second term. All other contributions make up the third, sub- subleading order term. 

The solution of Eq. (||) is straightforward, 

= &Aiexp(^-^^ [1 + 0(5)] , (9) 

where 6/i = 25 exp(&'/2). A more thorough analysis jsj shows that b is actually independent of S, and 
solely determined by b', i.e., by the contributions of screened gluons and the constants under the collinear 
logarithm. The constant d does not appear to leading and subleading order in the solution of the gap 
equation; it is part of the terms of order 0{g) on the right-hand side of (|9|). 

Observe that the power of the coupling g in the exponent is reduced as compared to the BCS result (^) 
or the result ^ for massive boson exchange. This was first noted by Barrois jo), but never made it into the 
published literature. Barrois assumed that the gluon propagator is the same as in vacuum and consequently, 
as discussed above, all gluon degrees of freedom are unscreened. 

In a dense medium, however, only almost static magnetic gluons are unscreened, while electric and non- 
static magnetic gluons are screened. In a seminal paper. Son derived the result (^ with the gluon 
propagator in the so-called "hard-dense-loop" (HDL) hmit 101^], and obtained 



Son also gave an estimate for the constant b, 



6=^, (11) 



with some constant 60 of order one. The parametric dependence of & on 5 arises from the gluon mass, Eq. 
(^. Several authors |^,||,|l2|,^ have confirmed the results ( [To| ) and (|l^) and refined the estimate for 60, 



60 = 512 TT* j b'^ , (12) 

where the dependence on Nf arises from that of the gluon mass (^, and 6q = 1 under the present approxi- 
mations. 
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In a dense medium, not only the gluon propagator is modified. Brown, Liu, and Ren included a finite, 
/x-dependent contribution to the quark wavefunction renormalization in their estimate for 6q, 



They also asserted that there are no further corrections to &q at this order in g, i.e., 6q = 1 + 0{g). 

Beane, Bedaque, and Savage [p^ argued that the coupling constant in the gap equation should not be 
taken at the scale /x, but at the scale of the momentum of the exchanged gluon. Using the standard running 
of the coupling constant in the vacuum they obtained 6q = exp[33(7r^ — 4)/64] ~ 20. At present, it remains 
an interesting open problem to check the validity of their arguments by systematically computing vertex 
corrections to the gap equation. Furthermore, their estimate could be improved replacing the vacuum 
running of the coupling by that in a dense medium. 

Sub-subleading contributions, i.e., of order g in the constant 6g, or equivalently, contributions to the 
constant d in Eq. (^, can arise from a variety of effects, for instance from the finite lifetime of quark 
excitations away from the Fermi surface |]l6|. Also, an apparent gauge dependence of the result (||) surfaces 
at this order However, the gap on the quasiparticle mass shell is a physical observable and thus in 

principle gauge- independent. This indicates that, for a complete solution to sub-subleading order, one has 
to go beyond the one-boson exchange approximation in the gap equation. This task appears, at least at 
present, too difficult to be feasible. 

In this paper, I focus instead on a contribution which is potentially of leading or subleading order. At 
small temperatures T ^ (/) ^ fi, the dominant contribution to the one-loop gluon self-energy comes from a 
quark loop; it is ~ g^^J■^, while gluon (and ghost) loops contribute a term ~ g^T^ and are thus suppressed 
p8| . In the standard HDL approximation, however, the quark excitations in the quark loop are taken to 
be those in a normal and not a superconducting medium. This is in principle inconsistent, as the fermion 
excitation spectrum in a superconductor differs from that in a normal conductor In Rcfs. p8| , p^ I 
argued that using the correct gluon self-energy in the gap equation could in principle affect the value of the 
color-superconducting gap . Consequently, I derived general expressions for the gluon self-energy in the 
two- and three-flavor case, and computed the self-energy in the limit of vanishing energy and momentum, 
which yields the Debye and Meissner masses in a color-superconductor. 

It turned out that in the two-flavor case, where the color-superconducting condensate breaks the SU(S)c 
color gauge symmetry to SU{2)c, the three gluons corresponding to the unbroken SU{2)c subgroup remain 
massless, like in the vacuum [nsi. Naively, if the propagator of these gluons is the same as in vacuum 
also at nonzero energy po and momentum p, this would change the number of modes mediating long-range 
interactions. Remembering the above discussion of the gap equation (^, this would affect the leading-order 
contribution to the gap equation and could in principle change the value ( |l0|) for the coefficient c. 

It turns out, however, that for zero gluon energy po = and gluon momenta |p| = p 3> 0, the self-energies 
approach their HDL form [p^ . Moreover, for energies and momenta much smaller than the gap, po, p <^ cj), 
these gluons do not behave exactly like in vacuum, but like in a medium with a large dielectric constant 
p]| . A change of the leading-order contribution to the gap equation, and thus of the coefficient c, appears 
therefore unlikely, but a change of the subleading contribution cannot be excluded. In order to decide 
this question, however, knowledge of the gluon self-energy in the limited range of energies and momenta 
POj P ^ 0; or for po = and (j) is insufficient; one has to compute the gluon self-energy for energies and 
momenta of relevance for the gap equation (^. As the integral on the right-hand side of (|l|) is dominated by 
quark energies and momenta close to the Fermi surface, the range of relevant gluon energies and momenta 
is Po, p fj,, but not the more restrictive range pq, p <^ (j) or the special case po = 0, p ^ (f>. In the present 
paper, I therefore compute the gluon self-energy for pq, p <C /i in a two-fiavor color superconductor. 

The gap equation in the two-flavor case was derived in Assuming that Cooper pairs are formed from 
quarks with fundamental colors 1 and 2, such that the condensate has fundamental color (anti-) 3, one 
obtains 



Here, is the gap matrix in Dirac space, V is the volume of the system, 7^ are Dirac matrices, A^^ is 
the gluon propagator, and S+ is the off-diagonal element of the quasiparticle propagator in Nambu-Gor'kov 
space. 




(13) 




Q 



(K-Q) S+(Q)7, 



(14) 
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Only gluons mediating between gapped quarks enter the gap equation, i.e., gluons with adjoint colors 1, 
2, 3, and 8. Gluons mediating between gapped and ungapped quarks, i.e., gluons with adjoint colors 4, 5, 6, 
and 7, do not appear. Moreover, the propagators for gluons of the unbroken SU{2)c subgroup are identical, 
A^^ = = ^33: but different from the propagator for the eighth gluon, Agg . 



I work in pure Coulomb gauge, 



where the gluon propagator for colors 1, 2, 3, and 8 assumes the form 



A':^,iP) = -Sab 

A°UP) = , 

^biP) = -Sab {5'= 



n°°(P) 



Here, n°°(P) is the 00-component and 



Ka{P)^7^{S,,-P^PJ)'^L{P) 



(15a) 

(15b) 
(15c) 

(16) 
I will 



the spatially transverse component of the gluon self-energy n^^(P), p = p/p. In the following 
(somewhat imprecisely) term n|J° the electric self-energy, and n|j^ the magnetic self-energy. 

Apparently, in order to settle the question whether the modification of the gluon self-energy in a color 
superconductor (as compared to the HDL approximation) changes the solution of the gap equation ([l^), 
one only has to consider the following four components of the gluon self-energy: H^]*, H^ S, H * ^ , and Hgg. 
The imaginary and real parts of these four components are computed in Sec. ||. In Sec. HI, I determine 
the spectral densities of electric and magnetic gluons as a function of gluon energy and momentum. Section 
[V discusses the effect of the modification of the gluon self-energy in a color superconductor on the solution 
of the gap equation. I conclude in Sec. ^ with a summary of the results. Unless mentioned otherwise, 
throughout this paper I work at zero temperature, T = 0. I use natural units and the metric tensor 
5^'' = diag (+,-,-,-). 



II. THE GLUON SELF-ENERGIES 



The starting point of the computation of the electric and magnetic self-energies for gluons of color 1 are 
Eqs. (99a) and (99c) of |l|l. For TV/ = 2 flavors, 

n;?(P) = 45^/y^ E (l + eie^kr-k^) 



1 1 \ ei £2 - Ci 6 



Po + ei+e2+i'n Po - ei - £2 + 2 6162 
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n'ii(P) = -2 5' / 7^ E (l-eie2ki.pk2.p) 

1 1 \ £1 £2 - 6 C: 



Po + £1 + £2 + Po — £1 — £2 + «7? / 2 6162 



(17a) 



(17b) 



Here, ki = k + p/2 and k2 = k — p/2, (pi = (j)^ is the gap function for quasiparticles (e^ = -1-1) or quasi- 

antiparticlcs (e^ — —1) with momentum k^, the variable 6 is defined as 6 = ^i^i — fJ-j and Ci = \/l^"+"W 
the quasiparticle excitation energy. The difference between Eqs. ( p^ and Eqs. (99) of is that the analytic 
continuation to real gluon energies po has been made explicit via the ii] prescription and that there are 
no terms ^ Ni = [exp(ei/T) -|- 1]^^, since iV^ = at T = 0. (Note that > 0i > in a superconductor.) 
According to Eq. (78c) of |18] , the self-energy for gluons of color 8 is 



u^^iP) = -ui;-{p) + -m''iP) , (18) 



^Hr(P) + i" 

where Hq'' is the standard HDL self-energy. Its components of relevance for the following are |p^ , |l8| 



5 



Po 



PQ + p-k + irj 



an r 



n*(P)^-m — l-(k.p 



47r 



Po 



Po + P • k + i?/ 



(19a) 
(19b) 



where mg is the gluon mass (|7|). According to Eqs. (102) of |18|, the electric and magnetic components of 
the self-energy 11'''' are 



n"V) = 4.9^7^ E (l + e,e.k,.k,) 



ei £2 - ^1 6 + 4>2 



n*(P) 



Po + ei + £2 + Po — - ^2 + ill J 
i'^'/(W ^ (l-eie2ki-pk2-p) 

ei,e2 = ± 



2ei £2 



1 



Po + ei + £2 + «?7 Po - ei - £2 + 



ei £2 - Ci 6 



2eie2 



(20a) 



(20b) 



^1 <(y2- 



The difference between Eqs. (17) and ( |20D is the relative sign in front of the term 

In the following, I compute real and imaginary parts of Eqs. ( p^ ) and ( pO| ) separately, using the identity 



—V i TT 5{x) 

X + ?>/ a; 



(21) 



for the energy denominators in Eqs. ( [l^ ) and (^0|), where V stands for the principal value prescription. I 
begin with the imaginary parts. 



A. Imaginary parts 

In the computation of the imaginary parts it is sufficient to consider gluon energies po > 0^ as 

Imn(— po,p) = — Imn(po,p) for all self-energies ( |l7|) , (^, (po|), which is easy to prove using Eq. (|2l|), 
and noting that 5{x) — 5{—x). It is instructive to first consider the imaginary parts of the HDL self-energies. 
For electric and magnetic gluons one obtains the standard result [p^ , ^ 

ImnO°(P)c.-^^m2^0(p-po) , (22a) 

3 2P0 /, Po 
- m — 1 - - 

4 " p \ p- 



Imn^(P)^-^^m2^ [l-^-^]e{p-p^) , (22b) 



corresponding to Landau damping for space-like gluons. 

For Imn5'5(P) and po > 0, the delta-function originating from the first energy denominator in Eq. (17a) 



has no support, since f-i > > 0. For the delta-function originating from the second energy denominator, 
use the fact that the interesting range of gluon energies in the gap equation (|lj) is po ^ ^, as the quarks 
in the Cooper pai r are close to the Fermi surface where ~ 0^ ^ ^. Then, the only term in the sums over 



ei, 62 in Eq. (17a) that one has to keep is the one where e-i — e2 ~ +1; for = —1, ~ |fci -I- ~ /Lt is far 
from the Fermi surface and the delta-function has no support for po ^ ^. 

Furthermore, the fermion momentum k in the loop is close to the Fermi surface, k fi, while the interesting 
range of gluon momenta in the gap equation is p <JC /e, which allows to approximate 

pk 



ki • k2 


~1 , 


^1 




6 





2 

P k 



(23a) 

e+ , (23b) 
e- , (23c) 



6 



where = k — fi. Finally, setting ~ (/)2 = (f), and defining e± = yCi + 0^: o^ic obtains 

Imn?;(P) ^-ng' [-^ 6{po - e+ - . (24) 

J (27r)-^ 2e+e_ 

First, note that since c-j- + e_ > 20, the delta-function has no support for pq < 2(f>: the imaginary part of 
the gluon self-energy vanishes below the threshold 20 for quasiparticle-quasihole excitations. This is clearly 
differ ent from the HDL self-energy which has an imaginary part for all gluon energies < po < p, cf. Eq. 

The integration over k is best done in spherical coordinates, choosing p = (0, 0,p). Then, the integration 
over the polar angle ip becomes trivial. Replacing the integration over azimuthal angle ^ by a; = cos^ 
and over fc by ^ = fc — /i, one can use the delta-function to either perform the integration over x, or over 
^. Choosing the latter, for a given value of x one finds that there are two roots of the argument of the 
delta-function, which differ by an overall sign. 



It is instructive to draw these solutions in the [x, i^)-plane, cf. Fig. |^. In part (a) of this figure, corresponding 
to Po > Ep, where 



Ep = v/p2+4 02 , (26) 

these two roots exist for all values of a; £ [—1, 1], while in part (b), corresponding to 2(p < po < Ep, they 
merge at ^ = and x = itpo/p, where 



; (27) 

Po 



there is no solution for 1 > > tpo/p. For pa < 20, there is no solution either, as already mentioned 
above. Note that Ep = 20 for p = 0, while Ep ~ p for p ^ 20. Analogously, t = for po = 20, while i ~ 1 
for Pq > 20. 

In the limit 0^0, the a:- integration decouples from the integration over ^, and by Eq. (p5[), the value of 
the latter is simply constant, = ±po/2. This is indicated by the dashed lines in Fig. |l. 
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tPo/2 

ULP 
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Po>E, 
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2(^<Po<Ep 



(a) 




(b) 
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FIG. 1. Solid lines: the function from Eq. 

the case = 0. 



for (a) po > Ep, and (b) 20 < po < Ep. Dashed lines represent 
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Evaluating the ^-integral at the roots (|25|), one obtains 



a2 /•" 



dy ■ 



PapJo (1- 2/2)3/2^^2 _y2 



(28) 



where w = min (t^p/po). The integral can be expressed in terms of elliptic integrals, yielding the final answer 



Po 



2 — p 
9{Ep ~~ Po) [E(t) - s2K(i)] + 0(po ~ Sp) 



Po 



pI~p^ 



, (29) 



where s = \/l — t'^ = 2 0/po, a = arcsin[p/(tpo)], and F{a,t) and E{a,t) are elliptic integrals of the first 
and second kind, 



F{a,t) 



dx 



1 



— sin^ X 



, £;(«, t) = / dx - 



sm X , 



(30) 



while K(i) = F{n/2,t) and E(t) = E{TT/2,t) are complete elliptic integrals of the first and second kind. At 
Po — Ep, IihHqq is continuous, since then a = tt/2 and the square root in Eq. ( |29| ) vanishes. 

It is important to note that the limit — > of the result ( |29| ) exists for all P: as (/> — > 0, a — > arcsin(p/po) 
and t ^ 1, such that E{a,t) — > p/po, while F{a,t) — > i ln[(po +p)/(po ~ p)]- Therefore, for po > Ep, 
the imaginary part vanishes like ^ (j)^ jp^. In the other region, po ^ Ep, when t — > 1, E(i) — > 1, while 
s2K(t) 0. In summary, the 0-^0 limit of Imn"5(P) is the imaginary part ( p2aD of the electric HDL 
self- energy, 



lim lmW{\{P) = luYW^{P) 

0— >0 



(31) 



The computation of the imaginary part of the magnetic self-energy for gluons of color 1 is completely 
analogous. With the approximation 



ki • p k2 • p ~ (k • p)2 , 



(32) 



valid when p ^ fc ^ /i, one derives 

Imn*i(P)~-7r^m2 0(po-2< 



4 

e{po - Ep) 



Po 
P 

pI 



0{Ep-po) 



Pi 



E(t) + s2^K(t) 



Po 



pI -p2 



(33) 



Again, Imll^]^ is continuous at po = Ep, and the limit — > exists for all P, 

lim Imn*ii(P) EE Imn^(P) . 



(34) 



The computation of the imaginary parts of the electric and magnetic self-energies (g^ is straightforward; 
the final result is 



lm&^\P)^~n^mle{po-2, 



Po 
P 

Po 



e{Ep - Po) E(i) + 9{po - Ep) 



E{a,t) 



P_ 

Po 



9 9 

Po -P 



Imn*(P) ~ - e{po '2 4>)- { 6{Ep - po) 



4 3 
e{po - Ep) 



pI 



l-^^{l + s') \ E(0-s2 1-2^^1 K(t) 



yo 



pI 



{l + s')\E{a,t)- 1 



pI\ P 



p- ) Po 



9 9 

Po-P 



1^2 



F{a,t) 



(35) 



(36) 
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As in the previous cases, these functions are continuous at po = Ep, and the Hmit — > exists; together 
with Eq. ( p^ one concludes 

hm Imn«0(P)^ImnS]0(P) , (37) 
hm Imn*8(P)^Imn*(F) . (38) 

Figure ^ shows the imaginary parts for a representative value of the gluon momentum p = 4 as functions 
of gluon energy po (in units of 2(j)). Comparing the self-energies in a color superconductor (solid lines) 
with those in a normal conductor (dashed lines) one observes the following two main differences. First, as is 
obvious from Eqs. (M), (ISSl), (ph, and (jsel), ImH™'* and Imn°°'* vanish below Po — 2(f), representing the fact 



that for such small gluon energies it is impossible to excite quasiparticle-quasihole pairs which would lead to 
non-vanishing self-energies. For the eighth gluon, however, the functions Imllgg'* do not vanish because of 
the admixture of the HDL self-energy, Eq. (|8|). Second, while the imaginary parts of the HDL self-energies 
vanish above the light cone po = the self-energies in a superconduc tor are never zero, although they 



rapidly fall off above po = Ep. As will be discussed in more detail in Sec. this indicates that, already at 
one-loop order, gluons are damped on the quasi-particle mass shell ||2^. Both features persist for all gluon 
momenta p. 

Another interesting property is that the electric self-energy for gluons of color 1 and the magnetic self- 
energy for gluons of color 8 are continuous, while the magnetic self-energy for gluons of color 1 has a 
discontinuity at pq — 2 cj) and the electric self-energy for gluons of color 8 has one at po = 2 and one at 
Po = p, the latter arising from the discontinuity of the HDL part in Eq. (p^). This leads to characteristic 
features in the real parts, as will be discussed below. Again, these features are generic for all gluon momenta 
P- 

I conclude the discussion noting that for po S> (/> the imaginary parts of the gluon self-energies approach 
the corresponding HDL expressions. They deviate from the latter only for po ~ cj). 



p=4(]) 
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Po/2(^ Po/2(^ Pj2<\> 

FIG. 2. The imaginary parts of the gluon self-energies as function of energy po for p = 4(j}. Figures (a), (b), and 
(c) show the electric self-energies, (d), (e), and (f) the magnetic self-energies. The solid lines in (a) and (d) are the 
imaginary parts of the self-energies for gluons of color 1, in (b) and (e) they are the imaginary parts of the functions 
n°° and ri' occurring in the self-energies for gluons of color 8, the imaginary parts of which are shown in (c) and (f). 
The dashed lines in all graphs are the imaginary parts of the HDL self-energies. 
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B. Real parts 



There are in principle two ways to compute the real parts of the gluon self-energies, either as a principal 
value integral from Eq. (|T]), or from the dispersion integral 



I r°° /I 

ReU{po,p) = -r dtjlmn(cj,p) 

Jo \UJ + 



po uj -po 



+ C , 



(39) 



where C is a (subtraction) constant, and where n(poiP) is assumed to be analytic in the upper complex 
Po-plane. 

In the HDL case, Eq. (|39| ) gives for electric gluons 



RenOO(po,p)^-3m^ 



fl-^ln 


Pq 


+ P 


) 


V 2p 


Po 


-P 





c: 



'00 



(40a) 



and for magnetic gluons 



Ren*(po,p) 



— m„ 
2 3 



pI 



pI 



Po 



In 



p^ \ p'- J 2p 
Comparison with a direct calculation via Eq. (|l|) deter mines 



/^OO _ 
Oq = 







Po+P 



Po~P 



(40b) 



(41) 



such that the term —2/3 in Eq. (40b) is cancelled by Cq and the magnetic self-energy has the correct zero- 
energy limit, RenQ(0,p) = 0, representing the absence of magnetic screening to one-loop order. Thus, the 
standard expressions for the real parts of the HDL self-energies are recovered JTo| , p3| . 

It should be remarked that the constant Cq is not a subtraction constant in the mathematical sense. A 
true subtraction constant would be required if ImnQ(po, p) did not vanish when pq oo. However, as can 
be seen from Eq. (|2^) , the imaginary parts of both electric and magnetic HDL self-energy are zero in this 
limit. The origin of a nonzero constant Cq is purely physical. As discussed in the last subsection, within 
the present approximations only quasiparticle-quasihole excitations near the Fermi surface give rise to the 
imaginary part of Hj^^. This holds also in the normal conducting phase, i.e., for IrnHp. When computing 
the real part via the dispersion integral ( ^9| ) these contributions are correctly taken into account. However, 
what is not accounted for is that, in contrast to the electric self-energies, the real parts of the magnetic 
self-energies also receive a contribution from quasiparticle-antiparticle excitations, cf. the discussion in Refs. 

There, I explicitly computed this contribution for po = P = in a color superconductor, and found 
it to be rrig, see also [ p5|j2^ ]. As antiparticles are always far from the Fermi surface, it actually does not 
matter whether one considers these excitations in a superconductor or in a normal conductor, or for nonzero 
gluon energy and momentum, provided po, p <^ ^, as is the case here. One therefore concludes that the 



constant Cq 



arises precisely from particle-antiparticle excitations. 



and that = 



0, because the 



electric self-energies do not receive such contributions. 

For the real parts of the gluon self-energies in a two-flavor color superconductor, I was not able to find 
expressions in closed form, so that the computation had to be done numerically. One can either directly 
use Eqs. (p7|), ( |2C| ) with Eq. (pi]), or the dispersion integral (^9|). It turns out that the second way is 
numerically far simpler; details can be found in Appendix This method does not specify the values of 



the subtraction constants and C* 
(p]), one immediately concludes 



1,8, but from the above discussion on the origin of the result 



00 _ r<oo — r<oo _ n 
1 = "-^88 = "-^o — 



(42) 

|2C|) with Eq. 



This has also been confirmed by a direct numerical calculation of the real parts via Eqs. (uTl 

®- n n 

Figure y shows the real parts of the gluon self-energies corresponding to the imaginary parts of Fig. g. 
The various discontinuities of the imaginary parts discussed above appear in the real parts as logarithmic 
singularities, due to the relationship between real and imaginary parts dictated by the dispersion integral 
(B3). Consequently, ReH™ has a logarithmic singularity at po — p, cf. Figs. (a), (b), (c), on account of 



the corresponding discontinuity of the imaginary part, while ReHg is regular, see Figs. ^ (d), (e), (f) 
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By the same line of arguments Re 11]']', shown in Fig. ^ (a), is regular. The peak at po ~ Ep visible in 
this figure is caused by the rapid decrease of the imaginary part in Fig. |^ (a), but it is not a true singularity 
as long as </> is nonzero. It becomes a singularity (the same as the one of Re 11°°) in the limit cf> 0. In 
contrast, Rell]] has a true logarithmic singularity at po — 2(1), Fig. |^ (d). This singularity persists for all 
nonzero (f), but obviously moves towards po ^ as cf) 0. 

As seen in Fig. ^ (b), the discontinuity of Imll^" ai pq — 2(f) causes a singularity in Re ft™. Together 
with the singularity of ReflQ" at po = P, this causes two logarithmic singularities in Rellg^, cf. Fig. ^ (c). 
As for Ren]], the peak at po = Ep is not a singularity, but is due to the sharp (but continuous) decrease 
of ImII°". Only in the limit 0^0, this becomes a true singularity and merges with the one of Rell™ at 
Po = p. The singularity at po = 2(j) persists for all nonzero (j). Finally, ReH* is regular, cf. Fig. ^ (e). From 
Eq. ( p^ ) and the regularity of Re Ilg one then concludes that the same holds for Ren^s- Fig- | (0- 

From the above discussion one concludes that 

lim Ren?°(P) EEERcn™(P) , (43) 
lim Ren|8(P)^Ren*(P) . (44) 

These equations hold for all P (even for the first equation, if, for po = p, one includes cxd as a possible value 
for limRen??). 

The functions Ren*] and Rellgg also converge towards the corresponding HDL expressions, except in a 
region of size ~ (j) near the logarithmic singularity at po = 2 (j), where the deviation from the HDL self-energies 
is large (infinite at po = 2 </)). Of course, in the limit — > 0, the point po — 2 (j) moves towards po — 0, and 
the size of the region where there are substantial deviations shrinks to zero. In this sense, the limit 0-^0 
for ReH]] and ReH™ still exists for all P and coincides with the corresponding HDL expressions. 

As for the imaginary parts, the real parts of the gluon self-energies in a color superconductor approach 
the corresponding HDL expressions when po S> (j). They deviate from the latter only for po ^ (j). 



p=4(t) 
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Po/2(^ Po/2(^ 
FIG. 3. Same as in Fig. H, but for the real parts of the gluon self-energies. 
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III. THE SPECTRAL DENSITIES 



The spectral densities are defined by the relations 



1 



TT 

When lnin"°'*(po, p) 7^ 0, the spectral densities are regular and Eq. m^) is identical to 



p°"(po,p) = - 



lninO"(po,p) 



[p2 _RenO"(po,p)]' + [Imnoo(po,p)]' ' 

1 Imn*(po,p) 

^ [pI - p2 _ Ren*(po, P)]' + [Inin*(po, P)]' 



(45) 

(46a) 
(46b) 



When Imn°'''*(po, p) = 0, for given gluon momentum p the spectral densities have simple poles, determined 
by the solution of 



[/-Ren°°(po,p)]^^_oo(p) 







in the electric case and 

[Po-^''-Ren*(po,p)]p,,_.(p)=0 
in the magnetic case, and are zero elsewhere, i.e., 

p™-*(po,p) = -^""-*(p){<^h-c."°'*(p)]-<^bo+^™^*(p)] } . 
The solutions Ci;™'*(p) to Eqs. (E^) define the quasiparticle dispersion relations, and 



are the residues at the poles pa = (jj°'''*(p). 



a(AOO^*)-i(po,p) 



dpo 



po=cd'"'-*(p) 



(47a) 



(47b) 



(48) 



(49) 



CM 




p/2(^ 



FIG. 4. Schematic plot of the (po,p)-plane. The HDL spectral density assumes the form (^6|) in regions la and Ha, 
and the form in regions lb, lib, and III. In a two-flavor color superconductor, the spectral density is always of 
the form (^) in regions Ila, lib, and III, and of the form (^^ in regions la and lb. 
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In order to understand the structure of the gluon spectral densities in a two-flavor color superconductor, 
it is instructive to first remember the structure of the HDL spectral densities in the (po^p) -plane, cf. Fig. ^. 
Below the light cone po = p (regions la and Ha in Fig. the imaginary parts of both electric and magnetic 
self-energies are nonzero, i.e., gluons are Landau-damped. Therefore, the spectral density is regular and 
). Above the light cone (regions lb, lib, and III), the HDL spectral density is of the form 
([48|), with the dispersion relations po — cjq°(p) for electric and pq = a;Q(p) for magnetic gluons. The 
explicit formulae for the electric {i.e. longitudinal) and magnetic {i.e. transverse) dispersion relations can 
be found in ||l0| , ^ ; they are shown graphically in Fig. ^ (dotted curves). In weak coupling, nig ^ 0, and 
since Wq°(0) = a;g(0) = irig, the electric and magnetic dispersion curves lie entirely in region HI of Fig. ^; 
in regions lb and lib the HDL spectral densities vanish. Figure || further clarifies the shape of the HDL 
spectral densities as functions of po for an exemplary value of p. 

On the other hand, in a two-flavor color superconductor, the imaginary parts of the self-energies never 
vanish for energies po > 2 4>, i.e., in regions Ila, lib, and III of Fig. ^ Thus, in these regions the spectral 
densities are always regular and of the form (^6|). Nevertheless, for pq > Ep (region III) the imaginary parts 
are small, ImH(po, p) '^-^ 0^- One can still define dispersion relations po = w(p) through Eqs. (p7|). These are 
shown in Fig. ^. In weak coupling, mg 3> (f>, these dispersion curves are nearly identical to the corresponding 
HDL dispersion curves, and consequently also lie in region III. (A peculiar feature is the dispersion curve for 
electric gluons of color 8, which has a negative slope for small momenta p ^ <f>. At present, I cannot exclude 
that this behavior is merely an artifact caused by neglecting the mixing of these gluons with the excitations 
of the condensate. This will be clarified in a future publication p7|.) 

Although the quasiparticles are not stable, like in the HDL case, they decay on a comparatively long time 
scale F"-*^ ~ [ImH(po, p)/po]~^ ~ The spectral densities resemble "smeared" delta-functions which 

peak at po = ^(p) and have a width ^ ImH(po, p) ^ (jy^ . This can also be seen in Fig. ^, where the electric 
and magnetic spectral densities in a color superconductor are shown as functions of po for a fixed value of 
the gluon momentum, p = mg/2. 



of the form (46 



magnetic 




FIG. 5. Dispersion relations for (a) electric and (b) magnetic gluons. 



: 8 <f). The different regions of Fig. M are 



indicated by thin solid lines. The dotted curves are the HDL dispersion relations. The solid curves are for gluons of 
color 1, the dashed curves for gluons of color 8. For magnetic gluons, the dispersion curve for gluons of color 8 and the 
HDL dispersion curve are visually indistinguishable. Note the additional dispersion branch below po — 2(f> (region la 
of Fig. mj for magnetic gluons of color 1. This branch merges with the continuum for po >2(j) above p ~ 1.1 mg. 
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0.5 1 1.5 0.5 1 1.5 2 



FIG. 6. Electric (a,c) and magnetic (b,d) spectral densities for gluons of color 1 (a,b) and color 8 (c,d), for a fixed 
value of the gluon momentum, p = mg/2, and rUg =8 0. The dotted curves represent the HDL spectral densities. 
In order to make the delta-function at the quasiparticle dispersion curves visible, I used the regular form (^6|) of the 
spectral density with a numerically small, but nonzero imaginary part. 



For pq < 2<p (regions la, lb in Fig. ^) the imaginary parts of the gluon self-energies in a color superconductor 
vanish and the spectral densities are of the form ( [48| ) . It turns out that only magnetic gluons of color 1 have 
a dispersion branch in this region, cf. Figs. I (b) and|(b). The origin of this branch, the properties of these 
gluons, and possible implications for deconfinement in the unbroken SU{2)c sector have been extensively 
discussed in pH . Here it is sufficient to note that for pg, p <^ (j), where the dispersion relation is approximately 
linear, the magnetic self-energy can be written as |2l[| 

2 

771 

^niP)^-7r^Pl^ (50) 



6( 



such that the magnetic gluon propagator becomes 



A\,{po,p)~— 2' (51) 



with the dielectricity constant 



The gluon dispersion relation is 



^ = l + 7d^- (52) 



^ii(p) - vp , (53) 

with the gluon velocity 

= ^ ~ \/6 — < 1 , (54) 

where the approximate equality and the inequality hold in weak coupling. The residue along the dispersion 
relation is 
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For the arguments presented in the next section, it is necessary to consider the hmit cj) ^ oi the spectral 
densities. In this Umit, regions la, lb, and lib shrink to zero. In region Ila (which becomes the complete 
region below the light cone po — p), the regularity of the self-energies IlJJ and n|g implies that the spectral 
densities ^^'^ Pss ^^'^ ^1^° regular. Hence, they converge smoothly towards the corresponding HDL 
expressions. 



lim 

0- 



npTAP)^Po'iP) (56a) 
hm p|8(P)-p*(P) , (56b) 

for P in region Ila. 

As discussed in the previous section, the real parts of the self-energies and H™ also converge towards 
the corresponding HDL expressions in the limit </>—!■ 0, except in a region of size ~ (j) near the logarithmic 
singularity at po = 2 (p. This singularity is rendered harmless in the spectral densities; the latter simply 
vanish at this point. The deviation from the HDL spectral densities is, however, still large near this point. 
The limit </> ^ nevertheless exists, since the point po = 2(p merges with the point po = 0, where the HDL 
spectral densities vanish, and the size of the region where there are deviations simply shrinks to zero. Hence 
also 



\im p{,iP) = pUP) , (56 
hm p°g(P)=pg°(F) , (56d) 

— >{) 



for P in region Ha. Finally, in region III the real parts of all self-energies converge smoothly towards the 
HDL expressions when (j) —^ 0. Hence, the dispersion curves in a color superconductor become identical 
with the HDL dispersion curves. Furthermore, the width of the "smeared" delta-functions in Fig. ^, being 
~ ImH ~ (/)^, goes to zero and simultaneously their height becomes infinite: they become true delta- functions, 
because 

5 /tt 

lim — — (57) 

is a representation of S{x). [To prove this, use Eq. (^l|).] Thus, Eqs. ( |5^ ) also hold in region HI, and hence 
in all regions of the {pQ,p)-plaiie which have nonzero measure in the limit — > 0. 



IV. THE EFFECT ON THE SOLUTION OF THE GAP EQUATION 

We are now in a position to assess the effect of the modification of the gluon self-energies in a color 
superconductor on the solution of the gap equation. Let us first perform the Dirac algebra in Eq. (p^); 
details can be found in |^. After projection onto positive and negative energy, as well as left- and right- 
handed particle states, and neglecting the contribution from the antiparticle gap, the result for the particle 
gap function is 



Q 



A[,iK^Q)--Al,{K-Q) 



9 

3-kq 



1 + kq 



1 + k • q {k- qf 



(58) 



Here, e+(0+) = •\/(q — pY + (0^)^i and the chirality index r, I of the gap function was omitted, since left- 
and right-handed gaps decouple and obey the same gap equation. To simplify the notation, in the following 
I shall also omit the energy index of the gap function, as well as the energy index and the argument of the 
excitation energy e+(^"'"). 
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Now approximate the terms involving k • q as 



1+kq 



1 



3-k-q l + k-q {k-qf 



-1 , 



(59) 



which is appropriate to leading and subleading order in weak coupling. In the resulting equation, add and 
subtract terms with HDL gluon propagators, 



6(f>{K) 



2 2Tv- 0(Q) 

3^ y 4^ 



3 0(Q) 

4 ^ T/ „2 _ 2 



AO"(if-Q)-A*(if-Q) 
A;?(ir-Q)-A!]0(i^-Q) 

A*i(if-g)-A*(if-g)" 



A™(if-Q)-A™(if-Q) 

A|«(if-g)-A*(if-g) 



(60) 



1 r 
9 



(61) 



For 5<j){K) = 0, we recover the standard gap equation with the gluon propagators in HDL approximation. 
To leading and subleading order in weak coupling, this equation is solved by the gap function derived by 
Son 0. At the Fermi surface, the value of the gap is given by Eq. (^) with c and h from Eqs. ( p^ ) and (p^), 
respectively. The additional term 5(t){K) is the correction to this solution. In the following, drawing on the 
results of the previous sections I determine whether this correction enters at leading order (modifying the 
value of c) , subleading order (modifying the value of b) , or sub-subleading order (no modification of c and h 
at all). 

First, rewrite the gluon propagators in Eq. ([l4|) in terms of spectral densities, see also Eq. (33) of [0, 



AOO(P) ^ -1 

1/T 



1/T 



dr eP"-" 



l-L/l POO 

A*(P)= / dref"^ / dujp\uj,p) \ [1 



nsiuj/T)] 



(62a) 
(62b) 



Here, I momentarily reverted to nonzero temperature until the Matsubara sum over fermionic energy = 
—i{2n + 1)ttT in Eq. ( |60| ) has been performed. Consequently, nsix) = [e^ — 1)^^ is the Bose-Einstein 
distribution function. The term — 1/p^ in the electric propagator cancels the contribution of A°"(P) when 
Pq ^ oo jl^ . This term is the same for the electric gluon propagators in a superconductor and the electric 
HDL propagator, because, as we have seen in the previous sections, there is no difference between these 
propagators for po ^ 0- 

To proceed, one also needs the spectral representation 



4>{Q) 



9o -e 



1/T 



2eo 



{ [1 - nFiejT)] e-^^^ - npiejT) e^'^ } , 



(63) 



where np{x) = (e^ + 1)~^ is the Fermi-Dirac distribution function. The Matsubara sum over go can now be 
performed in the standard manner. After analytic continuation onto the quasiparticle mass shell, — s- ek+irj, 
and abbreviating 0(eg, q) = (f>q the equation for S(j)k reads 



d^q 

(27r)3 2e 



q Jo 



SpZiu, p) - i Spf^iu, p) - 5p\^[u, p) + i 5pl^{uj, p) 



1 



1 



where p = k — q, and where I introduced 



"Paa — Paa Po 



00, t 



Eg - Efe 



a — 1, 



(64) 



(65) 



Neglecting a possible imaginary part of the gap function (cf. discussion in |3| ) , it is implicitly assumed that 
the energy denominators in Eq. ( |64[ ) are evaluated with the principal value prescription. 
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Finally, note that the spectral densities are isotropic, i.e., they only depend on |p| = p. Substituting 
the integration variable cos = k • q by p, and interchanging the order in which the p- and w-integrals are 
performed, one obtains 

<5^, ^ Aq ^ r dc. \jl,[u;) - J^',{^) ~ \ \Jk{^) " j ( ^ ' + — ^ ) ■ 

327r^ eq Jq y 9 ^ ' } \uj + Cq + uj + eq - J 

(66) 

where I defined 

rdpp6pZ'\u;,p). (67) 
Jo 

In deriving Eq. ( |6^ ) with Eq. ( |67| ) I also approximated k/q~ 1, fc + g — 2/i, |fc — q| — 0, which is valid, since 
the q-integral peaks at the Fermi surface, see also [|. Moreover, since the gap function (j)q depends strongly 
on g 0, it is permissible to restrict the momentum integration to a small region around the Fermi surface, 
/i — (5 < q < + 5 where 5 ^ fi. 

The important point is that the functions J^aa'* dimensionless, i.e., they can only be functions of 
the dimensionless variables Lo/nig and (j)/mg (or cj//i and (/>//i; additional factors of g turn out to be of no 
consequence for the following argument). Let us now estimate these functions. 



A. Estimate for J^f 



For u < 2(j), Pi? = 0, and Spil ~ —p'o'- When computing the p- integral in Jii, one samples in regions 
la and lb of Fig. 0. From the explicit form of the spectral density [p|jlO|j24|, 



p''o"{Lu,p) = 9ip-uj) 



2M^ uj 



TT p 



1 In 

2p 



P + OJ 



p — UJ 



2M^ - 

P 



(68) 



where M = "iumg/A, one deduces that Pq vanishes in region lb, while in region la peaks for p of order 
(and shghtly larger than) w, and then falls off ~ Neglecting powers of g, J^i is at most 



(69) 



For u! > 2ip, the p-integral samples regions III, lib, and Ila. In the latter, both pJi and Pq" a-re regular, 
consequently Spli is regular, and from the discussion in the previous section one concludes that it vanishes 
in the limit (/) — > 0. In principle, one could therefore Taylor-expand (Sp™ around (fi = 0. To leading order, 
one then has in region Ila 



r App5p\\{oj,p)^^ 



(70a) 



In region lib, Pq = 0, while p™ is regular and of order 1/m . The p-integral over this region is therefore at 



most 



, <ippSpll{uJ,p) 



(70b) 



Finally, in region III, (5pJJ is the difference between a "smeared" delta-function and a true delta-function, 
which is singular. However, the p-integral over JpJJ renders this singularity regular, and since Spqq vanishes 
in the limit (j) also in this region, the same argument as in region Ila applies. 



<^PP ^Poo(w,p) 



(70c) 
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Summarizing Eqs. (|69|) and ([70|), one concludes 

JiT(-) - -f (71) 

for all uj. 



B. Estimate for Ji-^ 



The estimate for the function J^^i is similar. The difference is that, in region la, there is also a contribution 
from yO*j. Explicitly, 



2^l 



dpp p*n(cj,p) = -pKi) Zl^ip) 



dp 



where the last term is the Jacobian from the argument of the delta-function in Eq. 
(^5|) one concludes 



/ dpp p[i{uj,p) ~ 1 . 



(72) 

With Eqs. (H) and 
(73) 



For our order-of-magnitude estimate, it is sufficient to use the approximate form of the magnetic HDL 
spectral density PJ2^] 



ujp 



valid for p ^ uj (this is where the magnetic HDL spectral density has its maximum). One then estimates 

r2M ,,, A 



(74) 



/ 

Jo 



dpp Pq{uj,p) 



Thus, to leading order 



Jliiuj) = / dpp Sp\^{uj,p) - 1 . 
Jo 



(75) 



(76) 



For uj > 2 (j), the situation is similar as for J'^^ and I refrain from repeating the details. One might worry 
that near uj = 2(1), where the deviation of p{i from pg is large due to the logarithmic singularity of ReH^^^j 
the above arguments do not hold. However, here Sp\i ~ — pg and with Eq. (|7j) one convinces oneself that 



still Ji*i 



/wg. In conclusion, to leading order 



Jl^{uj) ^ 9{2(j>-uj) X l + 9{uj-2(l)) X — 

m„ 



(77) 



C. Estimate for Jgf and Jgg 

The functions J^g and jTg'g are easy to estimate. For < 2 </>, ^Pgg' * is the difference between a spectral 
density with 2/3 of the value of the HDL self-energies and the HDL spectral density itself. Following the 
same line of arguments as in the previous cases, J7'8°g°'*(w) (j)/mg. For a; > 2 0, the situation is similar to 
that for the functions Cf^ '*. This immediately leads to the estimate 

^8°8"^*M-^. (78) 
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D. Estimate for 5(j}k 



We are now in the position to estimate d(j)k ■ With the exception of J^d in the region lo < 2 cj), all functions 
Jaa'* contribute at most terms of order (j)/mg to S(j)k- These contributions are suppressed by (jy/nig relative 
to the leading or subleading contributions. Additional powers of g are of no consequence, because they are 
accompanied by at least one power of (j) which is exponentially small in g, (/) ^ exp(— 1/g). Thus, these 
contributions cannot change the values of the constants c and b of Eqs. (^0|) and (p^). 

The only remaining contribution with a potential impact on the values for c or & is the one from . This 
contribution gives rise to a 6(l)k which is of order 

o(Pk^ g I M— In 

The integral can be estimated to be ~ 0, cf. Appendix ^ such that 6(j)k ^ 5^0- This means that 8(j)k 
contributes at most to the constant d in Eq. (||), i.e., to sub-subleading order to the solution of the gap 
equation. Thus, to leading and subleading order, the solution (H) with ( |lO| ) and ( |T^ ) is not modified. 

V. SUMMARY, CONCLUSIONS, AND OUTLOOK 

In this paper, I have computed the electric and magnetic self-energies for gluons of adjoint colors 1, 2, 3, 
and 8 in a color superconductor with two flavors of massless quarks. In an extension of a previous study 
which was restricted to gluon energies and momenta po = P = 0, as well as po = 0, ^ p ^ /i, 
I now considered the complete range poj p H- The results can be qualitatively summarized as follows: 
the presence of a color-superconducting condensate modifies the gluon self-energies as compared to the self- 
energies in a normal-conducting, dense medium only for energies po ^ (j). For po S> the expressions for 
the self-energies approach the standard hard-dense-loop results. 

This study was motivated by the following two observations. First, the weak-coupling result for the 
color-superconducting gap parameter depends sensitively on the number of unscreened (i.e. massless) gluon 
degrees of freedom. Second, in a two-flavor color superconductor not only magnetic but also electric gluons 
are massless in the static, homogeneous limit [Q. This led to the speculation as to whether the modification 
of the gluon propagator in a color superconductor could change the value of the gap parameter to subleading 
order in weak coupling, Eq. (||) with (^0|) and (|l2|). Using the results for the gluon self-energies, I estimated 
that this is not the case. The physical reason is that the change of gluon properties as compared to those in 
a normal conductor is limited to the rather small region po ^ 

There remains one unanswered question: the peculiar behavior of the dispersion relation of electric gluons 
with color 8 displayed in Fig. ^ (a). As discussed above, this could be due to the fact that the present 
calculation neglects the mixing of this gluon with the excitations of the condensate. This will be addressed 
in greater detail in a forthcoming publication . Since this phenomena is restricted to momenta p ~ </>, it 
is, however, unlikely that the conclusions regarding the stability of the value of the gap parameter to leading 
and subleading order will have to be revised. 
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APPENDIX A: NUMERICAL CALCULATION OF THE REAL PART OF THE GLUON 

SELF-ENERGY 



The real parts of the gluon self-energies can either be computed from Eqs. (p7|), ( [2C| ) with Eq. (|2l|), or 
from the dispersion integral (^). In both cases, this amounts to evaluating a double integral. In the first 
case, this integral runs over ^ = fc — /i and x = cos 9, while in the second case, one has to compute elliptic 
integrals numerically in addition to the integral over oj. 

The second way is, however, the simpler one. As the integrand falls off ~ '^/^'^ , one only needs to compute 
the integral up to some finite value Vl which is sufficiently large for the required numerical accuracy. One 
then divides the interval [0, SI] into N pieces of size 50, = Jl/iV, such that w„ — nSfl, n — 0, . . . ,N, and 
computes the integral over uj piecewise with a generalization of the mean- value theorem as 



1 ^ 

Ren(po,p) ^ - Vlmn(w;,p) in 

TT ^ ^ 



n=0 



(Al) 



with cij* g Wn-i-i] (for practical purposes, = {uin+i +w„)/2 is a convenient choice). In agreement with 
the principal- value prescription in (39), one has to make sure to avoid w„ = Pq for any n. The advantage of 
this method is that the imaginary parts can be tabulated prior to the computation of the real parts, such 
that the computation of the double integral is effectively reduced to the computation of two independent 
ordinary integrals. 



APPENDIX B: ESTIMATE OF THE INTEGRAL IN EQ. (|79|) 

In this Appendix, I estimate the value of the integral in Eq. ([79|), 



I(efe) = / dg ^ In 



(2 + e,)' 



f2 - f2 



(Bl) 



In order to see whether S(f>k contributes to the constants b, c, or d in Eq. (||), it is sufficient to perform this 
estimate at the Fermi surface, k = fi, i.e., tk — (p. Since the gap function is peaked at the Fermi surface. 



tlq < (t)tJ. 



J(0) < 



dq 



In 



-S £g 



(2. 



a21 



"^9 



In 



{2^ + e,f^^ 



2 1 



where = q — fi. Now substitute the integration variable |^] 

2/ = In — — , 



(B2) 



(B3) 



with the result 



1(0) <2(j> 



ln(25/0) 



dy In 



cosh^(?;/2) + 1 



sinh^(y/2) 

The right-hand side of this equation can be further evaluated, 

r-l 1 , „ /•ln(25/0) 



da; 1 + x /•"H^"/'?; 
1(0) < 2 <^ 2 / — In + / dy In 

I J4>/{2S) ^ Jo 



1 



cosh^(?//2) 



(B4) 



(B5) 



where x = e~y . An upper bound for the first integral is given by extending the lower boundary to zero. The 
resulting integral can be solved analytically and has the value 7r^/4. An upper bound for the second integral 
is given by extending the upper boundary to infinity. An integration by parts transforms the resulting 
integral into 
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dy y 



tanh(j//2) 
cosh^(y/2) + 1 



Substituting x = cosh(i//2), 



dx 



1 a;(l+a;2) 



In 



{x+\/: 



X^-l] < 4 



a;(l + a;2) 



ln(2a;) , 



(B6) 



(B7) 



as ^/x'^ — 1 < X ioT X G [l,oo). The last integral in (B7) can be solved analytically. Collecting the results, 
the integral 1(0) is bounded from above by 



I{(f>) <2(j) 



111 
12 



2 (In 2)2 



(B8) 



In conclusion, T{4>) is of order const, x (f>, and thus 50^ ^ 
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